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1.  INTRODUCTION 


For  a double  pendulum  subjected  to  a circulatory  force  applied 
at  its  free  end,  Herrmann  and  Bungay  [1]  determined  the  dependence 
of  the  critical  loads  of  divergence  and  flutter  upon  a tangency 
coefficient,  a,  which  specified  the  orientation  of  the  applied  force 
relative  to  the  instantaneous  deformed  configuration  of  the  system. 

It  was  assumed  in  reference  [1]  that  elastic  hinges  at  the  joints  in 
the  double  pendulum  exerted  linear  restoring  moments  c<j>^  and  c (<t>2_4>i  )• 
where  c denotes  the  stiffness  of  the  hinges  and  <|>^  and  $2  are  angles 
that  specify  the  configuration  of  the  system.  To  study  the  influence 
of  support  flexibility  on  the  state  of  stability  of  this  system, 
Suglyama,  et  al  [2],  assumed  that  the  linear  restoring  moments 
exerted  by  the  hinges  were  represented  as  c^<t>^  and  c2(<t>2-<t>i ) • They 
then  derived  expressions  for  the  critical  loads  as  functions  of  the 
ratio  of  stiffnesses,  Ci^. 

The  objective  of  the  present  investigation  is  to  consider  again 
the  effect  of  the  characteristics  of  the  support  on  the  stability  of 
the  system,  but  here  it  will  be  assumed  that  the  support  consists  of 
a platform  of  mass  M that  is  attached  to  a horizontal  elastic  spring 
of  stiffness  and  that  is  constrained  to  move  on  a smooth  horizontal 
surface.  As  a result  of  the  introduction  of  this  additional  mass  and 
degree  of  freedom,  the  system  to  be  studied  is,  of  necessity,  one  of 

1.  Herrmann,  G. , and  Bungay,  R.  W. , "On  the  Stability  of  Elastic  Systems 
Subjected  to  Nonconservative  Forces,"  Journal  of  Applied  Mechanics, 
Vol  31,  1964,  pp.  435-440. 

2.  Suglyama,  Y.,  Maeda,  S.,  and  Kawagoe,  H.,  "Destabilizing  Effect  of 
Elastic  Constraint  on  the  Stability  of  Nonconservative  Elastic 
Systems,"  Theoretical  and  Applied  Mechanics,  Vol  22,  University  of 
Tokyo  Press,  Tokyo,  1974,  pp.  33-45. 
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three  degrees  of  freedom.  The  frequency  equation  for  the  system  Is 
derived,  and  equations  for  the  critical  loads  of  flutter  and  divergence 
are  derived.  Stability  maps  in  the  load-tangency  plane  reveal  regions 
of  stability,  divergence,  and  flutter,  some  of  whose  boundary  curves 
are  strongly  dependent  upon  the  values  of  the  mass  and  stiffness  para- 
meters M and  k. 

2.  THE  EQUATIONS  OF  MOTION 

Consider  the  system  shown  in  Figure  1 that  comprises  a double 
pendulum  mounted  on  a platform  of  mass  M that  is  constrained  to  move 
on  a smooth  horizontal  surface.  The  configuration  of  the  system  is 
specified  by  the  coordinate  x.|(t),  which  locates  the  center  of  mass  of 
the  platform  relative  to  a fixed  origin  and  by  the  angles  $ (t), 
<t>2(t)  formed  between  the  vertical  and  each  of  the  two  bars  in  the 
double  pendulum.  The  double  pendulum  consists  of  two  rigid,  weightless 
bars  of  equal  length  £ and  carries  concentrated  masses  m^  and  m? 
located  at  distances  a^  and  a^  from  the  hinges  0 and  A,  respectively. 
The  hinged  joints  at  O and  A exert  linear  elastic  restoring  moments 
c<t>j  and  ),  where  c denotes  a torsional  spring  constant,  and  the 

horizontal  motion  of  the  platform  is  restrained  by  an  extensible  spring 
of  stiffness  k,  which  Is  fastened  to  a rigid  wall  at  C.  A force  of 
magnitude  P,  which  acts  at  an  angle  a^,  relative  to  the  vertical, 
where  a designates  the  tangency  coefficient.  Is  applied  at  the  free 
end  B of  the  pendulum.  This  applied  force  is  a circulatory  force 
when  a f 0 and  a conservative  force  only  when  a = 0.  Moreover,  it  is 
said  to  be  tangential  if  a = 1,  sub-tangential  If  0<a<l , super- 
tangential  if  a>l,  and  anti-tangential  if  a<0. 
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Under  the  assumption  that  the  displacement  and  the  angles 

4^  and  <J>2  remain  small  throughout  the  motion,  one  can  show  that  the 
expression  for  the  kinetic  energy  T for  this  system  is 


T = (M  + m1  + m^)  x|  + ~ (ny2  + n^Jt2)  <j>2  + (m^  + m?£)  + 

+ 2 m2a2^2  + V2V2  + m2£a2V2  O) 

The  potential  energy  V associated  with  the  extensional  spring  and  the 
elastic  joints  is 


V = |*xf  + \ c (2*f 


2<f>i  d>  2 + 


(2) 


Furthermore,  the  generalized  forces  associated  with  the  applied  force 
P are  easily  shown  to  be 

Qx1  = - Pot4>2 , Q-|  = - P£(a<f>2  - Q2  = - Pt(a-1  )$  . (3) 

Therefore,  using  Lagranges'  equations  in  the  form 


_d_  fil  \ _ 91  + 3V 
dt  lax^  dx-\  3X} 


± f3T  , 31  + 9V  _ 

dt  'W  ’ 3(()n  3<t>n  " V 


n = 1,2, 


with  equations  (1)  - (3),  one  obtains  the  equations  of  motion 


(m-ja^  + m2^)xi  + (m^  + m,,#,2)^  + (2c  - PH)^  + 

+ m fta2$  - (c  - aP£)4>2  = 0,  (5) 

m2d2Xl  + m2£a2$l  ' C*1  + m2a2^2  + “ Pl^  ' a^*2  = 0 • (6) 

It  is  desirable  to  express  equations  (4)  - (6)  in  dimensionless 

form.  This  may  be  accomplished  by  making  the  changes  of  dependent 

1/2 

and  independent  variables  x^(t)  = dx(t)  and  t = x!l(m/c)  , 

respectively,  and  by  defining 
kil2  Pit 

< = — » Q = — . M = mp,  mn  = myn  , an  = an£  , n = 1,2. 
Hence,  equations  (4)  - (6)  can  be  written  as 

(y-j  + y2  + y)x  + kx  + (y1a1  + y2)$1  + y2a2$2  + aQ  <p2  = 0,  (7) 

(y1a1  + y2)x  + (y-|0^  + y2)$-|  + (2  - Q)^  + y2«2$2  - (1  -otQ)4>2=0  (8) 

y2a2x  + y2a2$-j  - <j>1  + y2a^$2  + [1  - (1  - a)Q]<J>2  = 0,  (9) 

where  now  x = dx/dx,  etc. 

3.  THE  PURELY  DYNAMICAL  SYSTEM 

Suppose  next  that  the  masses  m^  and  m2  are  located  at  the  ends  of 
the  two  bars  in  the  double  pendulum,  i.e.,  at  the  points  A and  B, 
respectively,  in  Figure  1.  This  immediately  implies  that  a^=a2=l , 
so  that  equations  (7)  - (9)  become 

(y^  + y2  + y)x  + kx  + (y^  + l>2)^  + P2$2  + ctQ«f>2  = 0,  (10) 
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(yj  + y2)x  + (y1  + y^^  + .(2  - Q)<J>1  + y2$2  - (1  - aQ)^  = 0,  (11) 

y2x  + y2$.j  - <J>-|  + V 2^2  + ^ " a)Q^<t|2  = 0.  (12) 

These  differential  equations  with  constant  coefficients  admit  a 
solution  of  the  form 

x(t)  = X^e1^,  0»n(T)  = Xn+ieiwT,  n-1,2,  i = H)1/2,  (13) 


where  the  X 's  are  constants  and  oi  denotes  the  dimensionless  frequency 
J 

parameter  for  the  system.  Substitution  of  equation  (13)  into  equations 
(10)  - (12)  yields  the  set  of  homogeneous  algebraic  equations 


l [Amnu)2  - Cmn  - QDmn]Xn  = 0,  m = 1,2,3 
n=l 


(14) 


where 


Aii  = yl  +VJ2+T* » 


21 

Si 


Pl+Pgi 

v*2  * 


Cl!  * K- 
C21  ’ °' 
C31  * °’ 
°11  - 0, 

°2I  * °’ 
°31  ’ °- 
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*22 

S2 


y1+u2, 

yi+y2> 


= y 


2* 


C12  = °* 

C22  = 2’ 

c32  " -1 

0,2  ■ 0. 

°22  ’ -1’ 
°32  ■ °- 


A!3  = V 


’23 


” Uo  * 


A33  = P2’ 

C31  = °* 

C32  = _1» 

°33  = 1 * 

D1 3 = 

°23  = °* 

°33  = ■(1"a)* 


(15) 


It  is  now  easily  shown  that  the  frequency  equation  associated  with 
equation  (14)  is 
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where 


P0a>6  " p2(j4  + P4(Jj2  " p6  = °» 


(16) 


P2  = (2  + K)y1P2  + y(y1  + 5y2)-Q[y2(y1  + 2y)  + yp-j  (1  - a)], 

P4  = y^l  + k)  + y2 ( 1 + 5k)  + y - Q[(y1<  + 3y)(l  - a)  + y (3  - 2a)  + 

+ y2  (3  + 2k)]  + Q2[(y  + y.|)(l  - a)  + y2], 

pc  = k[1  - 3(1  - a)  Q + (1  - a)Q2].  (17) 

It  may  be  noted  that  the  value  of  the  determinant  of  the  dimensionless 
inertia  matrix,  A = (Amn),  is  det  (A)  = PQ  = y-^y.  1,e'’  ttie  Proc^uct 
of  the  three  fundamental  mass  parameters  associated  with  the  system. 

If  none  of  these  mass  parameters  is  zero,  then  the  inertia  matrix  is 
non-singular,  and  the  system  is  said  to  be  a purely  dynamic  system  [3]. 

For  the  sake  of  being  specific,  suppose  now  that  \i^  = 2 and 
y^  = 1 , so  that  equation  (17)  becomes 

PQ  = 2y,  p^  = 4 + 2k  + 7y  - 2Q[1  + y(2  - a)], 

p = 3 + 7k  + y - Q[(4  + 2k  + 3y)(l  - a)  + 5 + 2k]  + 

4 

+ Q2[(2  + y)(l  - a)  + 1], 

p6  = k[1  - 3(1  - a)Q  + (1  - a)Q2].  (18) 

3.  Ku,  A.  B.,  "On  the  Stability  of  a Linear  Nongyroscopic  Conservative 
System,"  Zeitschrift  fUr  angewandte  Mathematic  und  Physik,  Vol  20, 
1977,  pp.  986-991. 
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If  the  system  becomes  unstable  by  divergence,  the  critical 

loads  , j = 1,2,  can  be  computed  from  the  condition  of  vanishing 

frequency.  Putting  u = 0 in  equation  (16),  one  obtains  p,  = 0, 

o 

whence,  by  virtue  of  the  last  relationship  in  equation  (18), 

(1  - a)Q2  - 3(1  - a)Q  +1=0,  (19) 

b b 

provided  that  k f 0.  But  equation  (19)  is  identical  to  the  expression 

for  the  divergence  loads  of  a double  pendulum  attached  to  an  immovable 

base,  which  has  been  discussed  by  Herrmann  and  Bungay  [1].  Clearly, 

the  mass  M of  the  platform  on  which  the  double  pendulum  is  mounted  and 

the  constraining  extensional  spring  of  stiffness  k have  no  influence 

whatsoever  on  the  values  of  Q for  the  system  of  three  degrees  of 

bj 

freedom  under  consideration  here. 

The  flutter  loads  for  the  system  can  be  determined  from  the 
condition  of  the  merging  of  two  natural  frequencies  of  the  system. 
Because  the  frequency  equation  (16)  may  be  considered  a cubic 
polynomial  in  oj2,  the  condition  for  the  coalescence  of  two  of  its 
roots  is,  according  to  reference  [4], 

p0[4p0p4  * P2P4  * 27pOPI  ' 1SP0P2P4P6  * 4P2P6]  = °’  (20) 

If  pQ  / 0,  this  condition  becomes 

%p4  ' P2P4  + 2?P0P6  - 18p0p2p4p6  * 4p2P6  ’ °\  (21) 

1.  Herrmann,  G. , and  Bungay,  R.  W. , "On  the  Stability  of  Elastic  Systems 
Subjected  to  Nonconservative  Forces,"  Journal  of  Applied  Mechanics, 
Vol  31,  1964,  pp.  435-440. 

4.  Walter,  W.  W.  and  Anderson,  G.  L.,  "Stability  of  a System  of  Three 
Degrees  of  Freedom  Subjected  to  a Circulatory  Force,"  submitted  for 
publ ication. 


Substitution  of  equation  (18)  into  equation  (21)  yields  the 
following  sextic  for  the  critical  flutter  load  Qg: 

6 

l rO  = 0,  (22) 

n=0  n e 

where  the  rn's,  which  will  not  be  reproduced  here,  are  lengthy 
polynomials  in  the  mass  parameter  p,  the  stiffness  parameter  k,  and 
the  tangency  coefficient  a. 

4 . NUMERICAL  RESULTS  FOR  THE  PURELY  DYNAMICAL  SYSTEM 

Equations  (19)  and  (22)  may  be  solved  numerically  for  the 
critical  loads  of  divergence  and  flutter,  respectively,  for  given 
values  of  the  parameters  a,  p,  and  k.  The  objective  of  these 
numerical  calculations  consisted  of  the  preparation  of  stability 
maps,  i.e.,  plots  in  the  Qa-plane  for  selected  values  of  p and  k. 

The  curves  in  the  Qa=plane  obtained  from  equations  (19)  and  (22)  form 
the  boundaries  of  the  various  regions  of  stability,  divergence, 
flutter,  and  divergence-flutter.  In  order  to  identify  the  system's 
behavior  in  a given  region,  it  was  necessary  to  examine  several 
eigencurves,  i.e.,  plots  in  the  Qu)2-plane.  Three  typical  examples  of 
are  shown  in  Figure  2 for  (a)  a = 0.3,  (b)  a = 1.125,  and  (c)  a = 1.175 
with  p = 1,  k = 1/10.  In  these  eigencurves,  one  can  observe  the 
points  of  coalescence  of  various  modes,  which  indicate  either  the 
onset  or  cessation  of  a state  of  flutter,  and  the  intersection  of  the 
curves  with  the  load  axis,  which  indicates,  in  the  cases  depicted,  the 
transition  of  a frequency  from  an  imaginary  value  to  a real  value. 
Whenever  a frequency  for  a particular  mode  is  such  that  o.'?<0,  it  is 


Figure  2.  Eigencurves  for  p * 1,  k ■ 0.1,  and  (a)  a * 0.3,  (b)  a 
and  (c)  a = 1.175  (2  of  3). 


understood  that  this  mode  grows  exponentially  In  time  (divergence). 

If,  however,  it  happens  that  w2  has  a complex  value  for  some  mode,  the 
motion  in  this  mode  consists  of  an  oscillation  with  an  exponentially 
growing  amplitude  (flutter).  It  may  be  observed  from  Figure  2c  that 
u)2<0  and  uj2  oj2  are  complex  numbers  when  Q>5.615.  In  this  situation, 
the  various  modes  of  the  system  simultaneously  exhibit  instability  by 
both  flutter  and  divergence. 

Stability  maps  for  k = 1/10,  1,  and  10  with  y = 1 are  shown  in 
Figures  3,  4,  and  5,  respectively.  The  various  regions  of  stability 
and  instability  are  labeled  as  follows:  stable  region  (SR),  flutter 

region  (FR),  divergence  region  (DR),  and  flutter-divergence  region 
(FDR).  These  figures  reveal  that  the  shape  and  extent  of  the  primary 
flutter  region  is  strongly  dependent  upon  the  value  of  the  spring 
stiffness  coefficient  k.  For  relatively  weakly  supported  bases,  as  in 
the  cases  of  Figures  3 and  4,  the  flutter  region  is  fairly  large  and 
extends  beyond  the  range  of  the  tangency  coefficient,  a,  shown  here. 
However,  as  the  spring  stiffness  is  increased,  say  to  k = 10  as  in 
Figure  5,  the  flutter  region  contracts  significantly  in  size  and  is 
bounded  by  a closed  curve.  As  the  value  of  < is  raised  still  further, 
the  shape  and  location  of  the  flutter  region  in  the  Qa-plane  approach 
those  reported  by  Herrmann  and  Bungay  [1]  for  a double  pendulum 
mounted  on  a rigid  base. 

In  Figures  6 and  7,  the  stability  maps  have  been  plotted  for 
y = 1/10  and  y = 10,  respectively,  with  k held  fixed  at  k = 1.  Comparing 
Figures  6,  4,  and  7 one  sees  that  the  shape  of  the  flutter  region  is 
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Figure  6.  Stability  map  for  < = 1 and  p = 1/10. 


significantly  influenced  by  the  value  of  the  mass  parameter  p. 

5 . NUMERICAL  RESULTS  FOR  A QUASI-DYNAMIC  SYSTEM 

It  is  of  some  interest  to  examine  the  system  of  equations  (10)  - 
(12)  under  the  supposition  that  the  mass  of  the  base  is  negligibly 
small.  Since  the  determinant  of  the  dimensionless  inertia  matrix  A 
is  det  (A)  = p-jp2p,  if  is  obvious  that  det  (A)  = 0 when  p = 0.  There- 
fore, in  the  case  of  a massless  base,  the  system  under  consideration 
becomes  an  example  of  a quasi -dynamic  system  [3],  whose  motion  is 
described  by  the  differential  equations 

(p-j  + p2)x  + kx  + (p1  + p2)$-|  + P2$2  + clQ(J>2  = 0,  (23) 

(p^  + p2)x  + (p.|  + p2)^  + (2-Q)4>i  + p2$2  - (l-aQ)(fi2=  0(24) 

P2x  + p2^  - 4>-|  + P2$2  + [1  - ( 1 -ot)Q]d>2  = 0,  (25) 

which  are  obtained  from  equations  (10)  - (12)  upon  setting  p = 0. 
Moreover,  because  A is  singular,  it  follows,  according  to  reference 
[3],  that  an  "internal  constraint"  must  exist.  This  is  easily 
identified  simply  by  subtracting  equation  (24)  from  equation  (23),  the 
result  being 

kx  = (2-Q)<t>1  - <)>2  . (26) 

Using  equation  (26)  to  eliminate  the  x terms  in  equations  (23)  and  (25), 
one  finds  a pair  of  equations  of  the  form 
2 

l [Amn  $ + (Cmn  + QDmnH  ] = 0,  m = 1,2,  (27) 

i n n 


3.  Ku,  A.  B.,  "On  the  Stability  of  a Linear  Nongyroscopic  Conservative 
System,"  Zeitschrift  fiir  angewandte  Mathematic  und  Physik,  Vol  20, 
1977,  pp.  986-951. 


where  now 


An 

= (y-j  ♦ 

y2)(2  + k - 0), 

A = 
12 

(<  - 

l)y2  " 

*21 

= u2(2 

♦ ^ - q). 

A22  ' 

y2(tc 

- 1). 

C11 

= 2k, 

C12  = 

°n  ■ 

-tc. 

d.2 

C21 

= -tc. 

C22  = ’C> 

BZ'\  ' 

0, 

°22 

It  should  be  noted  here  that  the  new  dimensionless  inertia  matrix  A, 
some  of  whose  elements  Amn,  which  are  defined  in  equation  (28),  depend 
linearly  upon  Q,  is  not  symmetric.  In  addition,  the  value  of  its 
determinant  is 

det  (A)  = icp^p2(2  + tc  - Q)  . (29) 

Assuming  a solution  for  equation  (27)  of  the  type  given  in 
equation  (13),  one  can  verify  that  the  frequency  equation  for  the 
quasi-dynamic  system  is  simply 

P2u)4  - p4a>2  + p6  = 0,  (30) 

where  p2,  p^,  and  p^  are  obtained  from  equation  (17)  upon  setting 
y = 0.  Equivalently,  equation  (30)  can  be  expressed  as 

(A  + QHV*  + (D  + BQ  + iq2)u>2  + CQ2  + EQ  + F = 0,  (31) 


where 


A = vi^2  (2  + k),  H = - 0 = [(1  + <)p-j  + (1  + 5<)p2], 


t 


| 

.1 

r A 

- - Cv 

h 


B = (2  + <)  [(1  - a)p1  + u2l  + + (1  + k)p2.  I = (a  - 1)^  - \i^t 

C = k(1  - a),  E = -3tc(l  - a),  F = tc  . (32) 

It  may  be  remarked  at  this  point  that,  if  the  H and  I terms  were 
deleted  from  equation  (31),  the  resulting  frequency-load  relationship 
would  be  such  that  the  eigencurves  would  be  conic  sections  in  the 
Qu)2-plane,  as  reported  in  reference  [1].  However,  in  the  present 
investigation,  the  values  of  H and  I are  generally  different  from 
zero  and,  therefore,  the  eigencurves  are  not  conic  sections.  It  may  be 
anticipated,  then,  that  the  present  system  may  possess  some  flutter 
characteristics  that  are  unusual  relative  to  those  described  in 
reference  [1]. 

The  critical  divergence  loads  can  be  obtained  once  again  from  the 
condition  Pg  = 0,  i.e.,  from  equation  (19),  that  results  when  w = 0 
is  substituted  into  equation  (30).  On  the  other  hand,  for  the  flutter 
loads,  the  condition  for  the  coalescence  of  the  frequencies  obtained 
from  equation  (30)  is  p2  - 4p„p,  3 0,  whence 

4 c.  o 

I2Q-+  2( IB  - 2HC)Q3  + (B2  + 2ID  - 4HE  - 4AC)Q2  + 

+ 2(BD  - 2HF  - 2AE)Q  + D2  - 4AF  = 0,  (33) 

which  is  a quartic  polynomial  in  Q. 

In  the  case  of  p-j  = 2 and  p-|  = 1 , equation  (32)  becomes 


1.  Herrmann,  G. , and  Bungay,  R.  W. , "On  the  Stability  of  Elastic  Systems 
Subjected  to  Nonconservative  Forces,"  Journal  of  Applied  Mechanics, 
Vol  31,  1964,  pp.  435-440. 
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•A-fc. 


A = 2(2  + <),  H = -2,  0 = -(3  + 7tc) , 

B = 9 + 4k  - 2a  (2  + k),  I = 2a  - 3, 

C = k(1  - a),  E = -3k(1  - a),  F = k . 

Consequently,  equations  (31)  and  (33)  can  be  expressed  as 

2(2  + k - Q)u>-  - {3  + 7k  + [2  (2  + k)  - 9 - 4k]Q  + (3  - 2a)Q2}u)2  + 

+ k(1  - a)Q2  - 3k(1  - a)Q  + k = 0 (34) 

and 

(2a  - 3)  2Q4  - 2 [4 ( 2 + k)(1  - a)2  + 2(7  + <)(1  - a)  + 5 + 2k]Q3  + 

+ [4(2  + k)2(1  - a)2  + 4(13  + 6k)(1  - a)  + 31  + 34<  + 4k2]Q2  - 

- 2[2(<  + 2) (k  + 3) (1  - a)  + (2k  + 1)(7k  + 1 5) ]Q  + 9 + 26k  + 41k2  = 0, 
respectively.  (35) 

Solving  the  frequency  equation  (34)  for  co2  as  a function  of  Q for 
given  values  of  a and  k,  one  can  construct  eigencurves.  Two  such 
curves,  plotted  for  a = 1.05  and  a = 1 .25  with  k = 1 , are  shown  in 
Figure  8.  In  Figure  8a,  the  frequencies  of  the  two  modes  of  vibration 
coalesce  in  the  first  quadrant  of  the  Qw2-plane  at  Q = 1.01.  Conse- 
quently, the  system  becomes  unstable  initially  by  flutter.  As  the 
value  of  Q is  raised  even  higher,  flutter  eventually  gives  way  to 
divergence  which  persists  over  the  remainder  of  the  range  of  Q shown. 
In  the  case  of  a = 1.25  (Figure  8b),  the  onset  of  flutter  occurs  at 
Q = 1.02.  Flutter  persists  for  1.02<Q<2.78,  and  in  the  interval 
2.78<Q<4.01  the  system  is  divergent.  However,  for  4.01 <Q<5 . 02  the 
system  is  stable  once  more,  although  flutter  now  occurs  a second  time, 
in  this  instance  at  Q = 5.02,  as  the  result  of  a second  coalescence  of 
frequencies. 
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Figure  8.  Elgencurves  for  < * 1 and  (a)  a = 1.05,  (b)  a s 1.25  (2  of  2) 


Performing  computations  with  equation  (35)  for  the  critical 
flutter  loads  and  using  the  data  obtained  earlier  from  equation  (19) 
for  the  critical  divergence  loads,  one  can  plot  stability  maps  in  the 
Qa-plane.  In  Figures  9,  10,  and  11,  some  results  are  shown  for 
< = 1/10,  1,  and  2.  Although  these  figures  resemble  the  stability  maps 
in  Figures  3,  4,  and  5 for  the  purely  dynamic  system,  it  may  be 
observed  that  the  flutter-divergence  region  in  a stability  map  for  the 
purely  dynamical  system  becomes  a flutter  region  in  the  corresponding 
stability  map  for  the  quasi-dynamic  system.  In  addition,  the  divergence 
region  in  the  upper  right  corner  in  each  of  the  Figures  9-11  is 
replaced  by  a stable  region  in  the  quasi-dynamic  case. 


Figure  9.  Stability  map  for  k = 1/10. 
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